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probed is on the order of 1 mm 3 . 



BACKGROUND AND SUMMARY OF THE 
INVENTION 

1. Field of the invention 

The present invention relates to a method and an appara- 
tus to quantify the optical scattering and absorption 
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ABSTRACT 

We present a method and apparatus for local and super- 
ficial measurement of the optical properties of turbid 
media. The depth probed is on the order of 1 transport 
mean free path of the photon. The absorption coeffi- 
cient, reduced scattering coefficient and a phase func- 
tion parameter are computed from a single measurement 
of the spatially resolved reflectance close to the source 
(around one transport mean free path). Measurements 
on biological tissues can be achieved using a probe of 
diameter less than 2 mm, and the average volume 



properties of a turbid media. More precisely the present 
invention relates to a non-invasive measurement, over a 
small area of the sample surface. Local and superficial 
characterization of biological tissues in vivo is a major 
application of this invention. 

2. Related Background Art 

Different techniques have already been proposed to 
quantitatively determine the absorption and reduced 
scattering coefficients of turbid media 1 . Most of the 
n on -invasive methods are based on the measurement of 
spatially and/or temporally-resolved reflectance. The 
principle is as follows: the turbid medium is illuminated 
by a light source. The backscattered light is measured by 
one or several detectors. Different types of measure- 
ments are possible, depending on the time-dependence 
of the illuminating source: steady-state (continuous 
source), time-domain (short pulsed source) or frequency 
domain (amplitude modulated source). The present 
invention relates to the case of steady-state measure- 
ments, performed at different distance p between the 
source and the detectors. 

The range of p values is an important point to consider, 
when comparing different methods based on the reflec- 
tance. First, the probed volume of the turbid medium is 
related to the source-detector separation p. The larger 
the source-detector separation, the deeper the average 
depth probed. Second, depending on the range of p, dif- 
ferent mathematical processing must be used to obtain 



the optica] properties from the raw data. 
At least two cases must be distinguished. 
1 ) The first case corresponds to source-detector separa- 
tions larger than several transport mean free paths. For 
biological tissues 2 , this case correspond to source-detec- 
tor separation of typically larger than 2 mm. An analyti- 
cal form of the reflectance can be obtained from the 



method for the determination of the reduced scattering 
coefficient only, using an optical beam of oblique inci- 
dence. Moreover, their analysis does not include the 
effect of the phase function. Kessler et al. (Patents 
5,284,137 and 5,645,0619) proposed a method for local 
dye concentration and scattering parameters in animal 
and human tissues, based on spatial and spectral mea- 



dttrusion equation, ll trie absorption coefficient is suf- suremenis. However, iheir methods do not enable tu Llic 



ficiently lower than the scattering coefficient (typi- 
cally at least ten times). In such a case, the relevant 
optical properties are the refractive index, the absorption 
coefficient and the reduced scattering coefficient. The 
average depth of probing is on the same order than the 
source-detector separation. 

Such methods have been already published, and are the 
object of patents (Ref.l, Patent 517,987 Tsuchiya, 
Patent 5,676,142 Miwa et al.). 

2) The second case corresponds to source-detector sepa- 
rations close to one transport mean free path. For bio- 
logical tissues 2 , source-detectors separation range 
typically between 0.1 to 2 mm, and the average depth 
probed is on the order of 1 mm. Such small 
source-detector separations enables to measure locally 
the optical properties. 

In this case, light propagation can be modeled using 
Monte Carlo simulations. In contrast to the previous 
case 1 ), it was found that not only the first moment of 
the phase function must be considered, but also the sec- 
ond moment Precisely, the relevant optical properties 
are the refractive index, the absorption, the reduced scat- 
tering coefficient and a parameter y taking into account 
the two first moments of the phase function, and that we 
called the "phase function parameter*'. 
Different methods have been proposed for local charac- 
terization of turbid media, and in particular of biological 
tissues. Wang et al. (Patent 5,630,423) proposed a 



determine the absorption coefficient, reduced scattering 
coefficient and the parameter y. 

In the present invention we present a method and appa- 
ratus for the measurement of the absorption coefficient, 
the reduced scattering coefficient and the said phase 
function parameter y, from the spatially-resolved reflec- 
tance data at short source detector separation. The 
parameters determination procedure is based on the 
analysis of the reflectance, at short source-detector sepa- 
ration, we performed with Monte Carlo simulations. 
These parameters, which can be measured at different 
wavelengths, enable us to characterize turbid media, 
such as biological tissues. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig.l. Description of the principle of spatially-resolved 
reflectance measurement. 

Fig.2. Density probability function for the mean depth 
of scattering event. Case of detected photons at 
distances 0.3, 1, 1.5 mfp\ 

Fig.3. Examples of reflectance obtained with different 
phase functions. Case of matched refractive index 
(n=1.0). 

Fig.4. Monte Carlo simulations and fits on the form 
7?(p) = lA(p,ji; > Y) + R/B(^)] 2 . Case of 7^1.9, 
numerical aperture of the source and detectors = 0.37. 



A = 0.0647p 0324 exp(-0.161p), 
B = 0.18653 m - 0.8466 n a 2 + 1.836m 3 
Fig. 5. a. Basic description of the apparatus. Case of 
measurement with an optical fiber probe 

Fig. 5.b Basic description of the apparatus. Case of 
measurement with sources and detectors directly in 
contac t w ith the turbid medium, 



Fig.9. Plot of the parameter = — fR for different y values 
dp 

(1.0, 1.8,2.5), different reduced albedo a* (1, 0.95, 0.9) 
and for fixed p=l mm.. Mismatched refractive index 
n=1.4. 



Fig. 1 0 Plot of jR(a)' - jR{a'= 1) for different y values 
and reduced albedo a*, and for fixed p-1 mm. 
Mismatched refractive index n=1.4. 



Fig. 5.c. Basic description of the apparatus. Case of 
non-contact measurements with a ID or 2D detector, 
coupled to an imaging device. 

Fig. 5.d Basic description of the apparatus. Case of 
non-contact measurements with a scanning devices to 
illuminate and/or collect the backscattered light. 

Fig.6. Examples of the sample ending of the fiber 
optical probe. 

6a. simple arrangement for a single measurement. 

6b. arrangement for symmetrical measurements. 

6c .Arrangement for multiple measurements. 

6d. Arrangement for multiple measurement, using an 

multicore fiber. 

Fig.7. Example of a measured the spatially-resolved 
reflectance with the probe 2a. The measurement, 
performed on a microsphere suspension, is 
superimposed to a Monte Carlo simulation. The optical 
properties, computed from published water properties 
and Mie theory are: n=1.33, m=0.0004 mm" 1 , m'=1.0 
mm" 1 , Y=2.2. The calibration was performed with a 
siloxane sample of known optical properties. 

Fig. 8. Relation between the parameters R(p=l mm) and 
ldplnR(p— 1 mm)l and the optical coefficients and p.^ 
Case of 7=1.5 and 1.9. Probe of refractive index = 1.5, 
sample of refractive index = 1 A, optical fibers diameter 
= 200 urn, NA=0.37 (source and collection). 



DETAILED DESCRIPTION OF THE INVENTION 
1. Definitions 

The concept of spatially resolved reflectance is illus- 
trated in Fig.l. Consider light impinging on a turbid 
medium (through air or through a light guide), in a given 
solid angle co source . The spatially resolved reflectance 
R(p) is defined by the backscattered light power in a 
given solid angle o^ttcctop at a distance p from the 
source, per unit area and normalized by the source 
power. The distance p is referred as the source-detector 
separation (also in the case when light guides or imag- 
ing devices are used for the illumination and for the col- 
lection of the backscattered light). 

R(p) depends on the optical properties of the turbid 
medium, defined below. Note that R(p) depends also on 
the source and detectors characteristics, i.e numerical 
aperture and sizes- 
It is commonly admitted that the fundamental optical 
properties of a turbid medium is determined by the aver- 
age index of refraction n of the medium, the absorption 
coefficient n a , the scattering coefficient m, an< * the 
phase function p(6) where 0 is the scattering angle. The 
absorption coefficient |X a [m* 1 ] is defined as the probabil- 
ity of absorption per unit infinitesimal pathlength. The 
scattering coefficient flj [m~ l ] is the scattering probabil- 
ity per unit infinitesimal pathlength. The phase function 
p(0) is the density probability function for the scattering 



y i 



angle 0. The phase function is normalized as follows: 

1 = 2n f p(8)sin9<f9. (1) 
J o 

The rfi 1 order moment g n of the phase function is defined 
as: 



** = 2n \ P n (6)p(e)sin0Je (2) 
where P u (6) is the Legendnr polynomial of order n The. 



investigation, the present invention is directed to a novel 
approach where the volume probed is much smaller, on 
the order of 1 (mfp') 3 . This is achieved by using only 
small source detector separations, typically from 0.1 to 
2 mfp*. The lateral dimension of probing is limited to 
this range of distances. 

A model of photon migration in tissues was necessary to 



first moment of the phase function is also called the 
anisotropy factor, and is often simply noted g (=gi). It 
represents the mean cosine of the scattering angle. The 
reduced scattering coefficient Ji^' is defined as: 



Ms' = mO-g]) 



(4) 



The transport mean free path (or reduced mean free 
path) mfp* is defined as: 



The reduced albedo a' is the ratio: 



(5) 



(6) 



As a result of the present invention described in more 
details in the next section, it is also necessary to define 
another phase function parameter called y: 

Y=(l-g 2 )/a-gl) (7) 
All the parameters listed above are referred as optical 
properties. They are wavelength dependent, and can 
vary in space and time. 

2. Reflectance measurements at distances close to 
one transport mean free path 

The methods described in the Patent 517,987 (Tsuchiya) 
are based on measurements with large range of 
source-detector separations, typically from 1 mfp* to 10 
mfp*. In such cases, the volume probed is on the order 
of 10-1000 (mfp 1 ) 3 . In contrast with such a large scale 



predict the relationship between the measured reflec- 
tance and the optical properties. Analytical solutions 
from the diffusion equation are not appropriate in our 
case because we are interested in the reflectance close to 
the source, at distance comparable to the transport mean 
free path [mfp*]. This is part of the invention to have 
performed Monte Carlo simulations to predict the mea- 
sured reflectance of an homogeneous semi-infinite tur- 
bid media. 

The exact diameter of the illuminating and collecting 
fibers, as well as their numerical apertures, have been 
taken into account in the simulations. The mismatch of 
index of refraction at the surface of the medium have 
been also taken into account in the simulation, by using 
the Fresnel law for each photon reaching the surface. 
This is also one results of this simulation to compute the 
average depth of probing, illustrated in Fig.2. It is dem- 
onstrated that only the superficial part of the turbid 
medium is probed if small source-detector separation 
are used. The average scattering depth of scattering 
event was recorded for each photon detected in Monte 
Carlo simulations. We present in Fig.2 the probability 
density function of this quantity. Fig.2 shows that the 
average depth of scattering is approximately around 

1 mfp\ Moreover it showed that the part located below 

2 mfp* were not playing a significant role in the mea- 
sured signal (for p < 1 .5 mfp*). 

The spatially resolved reflectance R(p), with short 
source-detector separations is more complex than in the 



ni 



case of large source-detector separations. Indeed, for 
small source detector-separation conditions, the inverse 
problem, i.e. calculating localized absorption and the 
reduced scattering coefficients, is necessarily sensitive 
to the scattering phase function. It is part of the present 
invention to have shown, from Monte Carlo simulations, 
that only the first and second moments of the phase 



extracted from measurements at short source-detector 
separation: ^i a , m' and y. 

It is important to note that the parameter y, cannot be 
estimated if large source detector separations are used. 
Indeed, the reflectance is sensitive to y, only at short 
source-detector separation. 

It was also derived from Monte Carlo simulations that 



function must be taken into account. Moreover, it was 
established that the influence of the these two moments 
are not independent. Indeed, it is a merit of this inven- 
tion to demonstrate that only one new parameter y, 
called "phase function parameter", which depends on 
the first and second moments can be used to characterize 
correctly the reflectance of tissues at short distance. Pre- 
cisely, y is defined as follows: 

y=(i-g 2 ya-gi) (» 

More, precisely, Fig.3 illustrates the fact that the param- 
eter y is the only predominant parameters of the phase 
function that must be taken into account (and not gj as 
frequently mentioned in literature). Reflectance curves, 
obtained from Monte Carlo simulations are shown in 
Fig.3. Four different phase functions were used for the 
simulations. Three phase functions are characterized by 
y^l.25. For comparison, the results for a fourth phase 
function with y=2.25 is also presented. For distances 
PH S * <0.3, the parameter y appear clearly to be the 
important parameter of the phase function. 
The parameter y depend on the characteristics of the set 
of scatterers, inside the turbid medium: shape, distribu- 
tion of sizes and refractive indices (medium and scatter- 
ers). Thus, the determination of y can provide important 
information about the investigated sample. 
In most cases, the refractive indices n and the source and 
detectors characteristics are fixed and known. In such 
case, a reduced set of three parameters only can be 



the r eflectance R(p), foi p ^lmfp r , can be r easonab ly — 
well approximated by the expression: 

R{p) = [A(p, u.;, y) + H/B^)] 2 (9) 

it is important for the description of the invention to evi- 
dence the fact that the function A depends on p, the scat- 
tering properties (i.e. \x s * and y) but not on m. In 
contrast, the function B depends on p^ but neither on the 
phase function (i.e. y) and nor p. An example of Equ.(9) 
is shown in Fig.4. 

This formulation is of great help to solve the inverse 
problem, as described in section 5.3. 

3. Apparatus 

In the first embodiment, The apparatus can be divided in 
three parts, described in Fig.5.a. 

The first part is the illuminating system. Different fight 
sources can be used: a) white sources, such as halogen 
or xenon lights, metahalides or fluorescent or phospho- 
rescent sources, b) the sources described in point a) 
where monochromators, filters or interference filters are 
added to select a given set of wavelengths c) sources 
such as laser, laser diodes, light emitting diodes or 
superluminescent diodes. 

In the first preferred embodiment, The light power is 
conducted to the investigated sample by the probe, 
which is the second part of the apparatus. The probe is 
preferably made of optical fibers, to illuminate and to 
collect the backscattered light. Different possible 



# 



arrangements are illustrated in Fig. 6. Two different 
modes of measurements can be chosen. First, one fiber 
is used to illuminate the sample and at least two others 
are used to collect the backscattered light at two differ- 
ent distances. Second, one fiber is used to collect the 
light and at least two others fibers, at two different dis- 
tances from the first one, are used to illuminate sequen- 



and fibers, light pipes or grin rods are also included in 
the present embodiment 

The third embodiment is described in Fig. 5.c.and Fig. 
5.d. The only difference with the first and second 
embodiment is that no contact measurements are per- 
formed. A collimating system allow for point illumina- 
tion on the turbid medium. An imaging system enables 



dally the sample. 

The arrangement of the different fibers can be replaced 
by any imaging system or image guide, such as multi- 
core optical fibers. 

The light collected by the probe is analyzed by the 
detection unit, which is the third part of the apparatus. If 
wide spectral light sources are used (such as halogen or 
xenon lights), a spectrograph should be put between the 
probe and the detector to get wavelength dependence of 
the backscattered signal (either in the source or detec- 
tion unit). Different types of detectors, such as photo- 
diodes, avalanche photodiodes or photomultipliers can 
be assigned to each collecting fibers. Simultaneous 
detection of each collecting fiber can also be achieved 
using linear or two-dimensional detectors such as 
Charge-Coupled Detectors (ID or 2D) or array of photo- 
diodes. 

The second embodiment is described in Fig. 5.b. The 
only difference with the first embodiment is that option- 
ally no optical probe, based on the use of fibers or light 
pipes or grin rods is used. The light source unit is 
directly in contact with the turbid medium, as well as the 
detector unit. Collimating optics, microoptics or imag- 
ing optics (DOE for example) can be put between the 
turbid medium and the actual light sources and detec- 
tors. The different type of sources and detectors cited in 
example for the first embodiment can be used for the 
second embodiment. Hybrid design, such as arrange- 
ments involving both direct contacts sensors or detectors 



to measure the spatial OistnbuUon 01 trie reflectance. 
The detectors can be either an array (ID or 2D) of detec- 
tors (Fig. 5.c),or a single detector (Fig. 5.d). In this last 
case, an scanning device is used to obtain the spa- 
tially-resolved reflectance. A fiber bundle, mulucore 
fiber or relay optics (grin rod or multiple lenses) can be 
put between the focal point of the imaging system and 
the detector(s). The different type of sources and detec- 
tors cited in example for the first embodiment can be 
used for the third embodiment 

It must be also noted that multiple measurements at dif- 
ferent locations allows to obtain an image of the differ- 
ent parameters u a , m* and y (the resolution is on the 
order on the mean source-detector separation used). 

4, Normalization and calibration 

The differences of transmission between each fiber are 
corrected using a measurement on a turbid phantom illu- 
minated uniformly. 

The background light, measured with the light source 
turned off, must be subtracted from the signal. 
In order to perform absolute intensity measurements, 
calibration is performed on turbid medium of known 
optical properties. Examples of such media are: 1) solid 
turbid medium which properties have been measured by 
other standard techniques 2) water suspension of micro- 
sphere of know size distribution and refractive index. In 
case 2) the scattering properties are calculated from Mie 



theory, and the absorption coefficient is assumed to be 
equal to the water absorption coefficient. 
A Monte Carlo simulation is performed with the optical 
properties of the calibration sample. The experimental 
reflectance performed on the calibration sample is mul- 
tiplied by a factor determined so as to fit the experimen- 
tal to the simulated reflectance. This factor is defined as 



The slope of the square root of R(p) is given by: 

( m i \ 



(ID 



5.3. Inverse problem 
Wc developed different methods to solve the inverse 
problem, which consists in extracting optical coeffi- 



the calibration factor. Each new measurement is multi- 
plied by the calibration factor. 

5. Signal processing 

5. J . Control of the homogeneity of the area probed 
Artifacts during a measurement, for example due to bad 
contact between the probe and the sample, or heteroge- 
neity of the sample, can be detected by the following 
procedure. Two illuminating fibers are disposed sym- 
metrically in regard to the collecting fibers (see Fig .6b). 
If the sample is homogeneous, the reflectance curve 
should be identical with either illuminating fiber. There- 
fore, comparing the two curves tests the heterogeneity 
of the investigated tissue region or detects obstructions, 
beneath the fibers. If the two curves are sufficiently 
close, the measurement is validated and the average of 
the two curves is calculated. 

5.2 Smoothing procedure and computation of the 
derivative of the curve. 

Functions in the form m x p exp(m 3 p) were found to fit 

always well Monte carlo simulation results for restricted 

range of distances. Smoothing of the experimental 

reflectance R(p) is obtained by fitting R(p) to 

m-ip exp(m 3 p). 

The determination of the slope of the logarithm 
(ln/?(p, p 5 \ a', y)) is also derived from this fit, using 



5p 

the following formula: 



^(ln*(p,|i,\a\Y)) = y + m 2 . 



(10) 



cients from the reflectance data. 
Method \. 

Monte Carlo simulations show that measurements of the 
reflectance intensity R(p) and the slope of lnR(p), deter- 
mined at a fixed distance p, can be used to derive u. s ' and 
ji a for a given y value. Fig.8 shows graphically the rela- 
tionship between fig' and ^ and the two parameters 
R(p=l mm) and ia p lnR(p=l mm) I for y= 1.5 and 1.9. 
We see in Fig.8 that m' ^ can not bc determined 
uniquely if y is unknown. Further insight for optimizing 
the inversion strategy is provided by three additional 
features of Fig.8. 

First, the determination of is only weakly influenced 
by y. Indeed in Fig.8 the errors induced by error in 7 are 
typically 10% for m' around 1 mm" 1 . Second, although 
absolute determination of m is not possible when y is 
not precisely known, relative absorption changes can be 
still evaluated. Third, the indetermination of the parame- 
ter y may be resolved by the values of R(p) and/or 
l3 p lnR(p)l at other distances. Therefore the following 
procedure can be used: 

(1) deterrnination of and M^a frorn R(P == 1 mm) and 
!3plnR(p=l mm)l for a set of values y. For example: 
7^1.0, 1.1, i.2,...,2.5. 

(^simulations with the different sets of and m 
obtained 

(3)comparison between the simulations and the reflec- 
tance profile for distances 0.3<p<2 mm. 
This last step allows us to determine the correct values 



of Y, and p^. Points 1 to 3 can be done iteratively 10 
evaluate y more precisely, using a finer discrimination of 
y values. 

Method 2. 

The inverse problem can be optimized, considering 
properties of R(p) of Equ.(9). 



value of R(p), which highly depends on u^. For given 
Hs\ y and p values, the dependence of reflectance on ^ 
is obtained by Monte Carlo results. From Equ.(9) we 
have: 



(14) 



where / 



Indeed, it can be^denveq rrom bqu.(^) triat trie quantity 
~^jR{p r u./, V> a > does not depend on the absorption 
coefficient u^: 



^(p,M/, Jlfl ,Y)=2^(p ) n;,7) (12) 
This property is confirmed in Fig. 9, where the quantity 



2^J R (P> M- 5 '> Y), evaluated at p=l mm, is plotted as a 
function of for y=l, 1.9 and 2.5 and reduced albedo 
0.95 and 0.9. In Fig.9, the parameter 



a'=l, 

a 



/R(p) , calculated from the experimental reflec- 



^-jR(p t u./,y) clearly depends on m' and y. In con- 
trast, the dependence on a* is almost negligible. 
Therefore, y and can be derived from the parameter 

dp 

tance R(p) (see section 5.2). Simultaneous determina- 
tion of y and m* require values of ~^jR{p) at different 
distances p (at least two). If y or is already known, 
the determination of the unknown parameter can be 

obtained from the value of ^-jR(p) at a single distance. 

a P B 

Convenient analytical approximation of ^-jR(p) can 

be obtained by fitting Monte Carlo results to polynomial 

functions. For example, in the case of fixed y= 1 .9, we 

obtained: 

Hs*=0.9 1 62-52.89x+ 1 795x 2 -l 8155x 3 +65428x 4 (13) 

where x= -~-+/R(p) at p=l mm, expressed in [mm -2 ].. 
This results is valid for a probe of refractive index of 1 .5 
and a sample of refractive index of 1.4, optical fibers of 
diameter 200 |xm, NA=0.37 (source and collection). 
Once fig* and y are calculated from the procedure 
explained above, is calculated from the absolute 



ulation. Particularly, they can fee well approximated by 
polynomial functions. For example, for y=1.9, probe 
refractive index of 1.5, sample refractive index of 1.4, 
optical fibers of diameter 200 u:m, NA=0.37 (source and 
collection): 



(15) 
(16) 



/ =-0.002257-8. 17 lm , +268.8u^' 2 
h = 0.0131 l+0.05184x-0.01974x 2 + 
0.003217x 3 -0.0001992x 4 

where x = ^W'/ j (unitless). 



Method 3 

Equ.(9) also show that relative measurements of u^, i.e. 
variation of |X a from a known value, is possible by the 
monitoring variation of the parameters jR(p) 
Indeed, we have: 



V*(p,n;,y, u fl ) - V*<P> Y. H fl = 0) (17) 

= R;5(n fl )-u;B(u a = o) 

This relation is illustrated in Fig. 10 in the case of p=l, 
and for the same three y values used in Fig.9 and for 
a*=l to 0.83. Fig. 10 confirms that the influence of y is 
weak in the quantity 



JR(P* Vs> 7. M - J R (P> Vs> T. M- fl = 0) . For known m\ 
the function B(a*) allows for the determination of a rela- 
tive absorption change Au. a =|i a - ii^, from a known 
value u-^j. Fig. 10 illustrates the case fi ao =0, but any 



hi 



other value of is possible. The interesting point is 
that B(a') does not depend on the phase function. 

We claim: 

1.A method for local and superficial (on the order of one 
transport mean free path) characterization a turbid 
media using all the following parameters: 



2 i 

The expression H = m : p exp(m 3 p) gives a smooth 
function R(p). 

The slopes ^Lv7?(p) and A(lnrt(p)) , or any mathe- 
matical combinations of these two latter quantities and 
R(p), can be obtained directly from analytical functions 
using the parameters m 2 or by numerical proce- 
dures from the function R = m x p exp(m 3 p) . 



1) the refractive index n of the turbid media 

2) the absorption coefficient |i a of the turbid media 

3) the reduced scattering coefficient p. s ' of the turbid 
media 

4) the phase function parameter y of the turbid media 
and comprising the step of: 

-measuring the spatially-resolved reflectance R(p) or 
any quantity which allow for the indirect determination 
of the said spatially-resolved reflectance R(p). 
-mathematical processing to compute at least one of the 
said parameter (n, y) or the relative variation of 

these said parameters. 

2.The method of claims 1, wherein said spatially 
resolved reflectance is measured by an optical fiber 
probe as described by the first embodiment in section 3. 
3The method of claims 1, wherein said spatially 
resolved reflectance is measured by a device put directly 
in contact to the turbid media, as described by the sec- 
ond embodiment in section 3. 

4. The method of claims I, wherein said spatially 
resolved reflectance is measured by non contact detec- 
tion unit, as described by the third embodiment in sec- 
tion 3. 

5. The method of claims 1 to 4, wherein the following 
processing is performed: 

- fit of the measured reflectance R(p) to the function: 

mjp 2 exp(m 3 p) (18) 
This fit give values for the parameters m lr m 2 and nvj. 



6. The method ol claims i to 5, wherein the absorption 
coefficient the reduced scattering coefficient fV and 
the phase function parameter y are determined on the 
basis of a fit of the spatially-resolved reflectance R(y, 
|i s \ u^, p) to a Monte Carlo simulations, or functions 
approximating Monte Carlo simulations. 

7. The method of claims 1 to 5, wherein the absorption 
coefficient the reduced scattering coefficient \i$ and 
the phase function parameter y are determined on the 
basis of the following form of the reflectance; 

R(p) = (A(p, y, u, s ') + m'B(m)) 2 (19) 

where the function A(p, y, and B(\i^) taking into 
account the sources and detectors characteristics and the 
refractive index of the sample. 

8. The method of claims 1 to 5, wherein the reduced 
scattering coefficient u^' and the phase function parame- 
ter y arc determined on the basis of the quantity 
S- JR(P> H/» y) » which depends weakly on the absorp- 
tion coefficient p^for 0.3<pfi s , <5. 

9. The method of claims 1 to 5, wherein the absorption 
coefficient m is determined based on the equation: 

rjRte)-/(y,u $ ')1 



(20) 



where / and h are functions that approximate Monte 
Carlo simulations. 

10. The method of claims 1 to 5, wherein variation of 
the absorption coefficient m is possible computed from 
the variation of the parameters jR{p) . 

11. The method of claims 1 to 5, wherein the m value 



I 



can be obtained from a function B, derived from Monte 
Carlo simulation, such as: 

jR{p, u;,Y.H 0 )-7*(P. ^ V, H fl *) (21) 

where jR{p, \l s \ y, u fl *) is a obtained from measurement 
on a sample with a known |i a * value. 
T2. The method of claims 1 to 6, wherein said turbid 
medium is a biological medium. 
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The spatially resolved reflectance of turbid media is studied at short source-detector separation (around 
one transport mean free path), using Monte Carlo simulations. The role of the phase function is carefully 
assessed. Particularly, the importance of the first and second moment of the phase function is 
demonstrated for distances between 0.5 to 5 transport mean free paths, whereas the effect of moments of 
higher order is negligible. Similarity relations are tested and are found efficient to reduce the number of 
relevant parameters, Indeed, only the following parameters must be t a ken i n t o an/mu it: i I i k i r . rim iivr 



index, the absorption coefficient, the reduced scattering coefficient, and a phase function parameter y 
depending on the first and second moment of the phase function. Approximate analytical form of the 
reflectance at short source-detector separation are also derived from the simulations. It is based on the 
follwing findings for source-detectors separation between 0.5 to 2 transport mean free path: first the 
slope of the square of the reflectance depends weakly on the absorption, and second, differences of the 
square root of the reflectance weakly depends on the phase function parameter y. These results give clues 
for the determination of scattering and absorption properties, using short short distances measurements. 



1. INTRODUCTION 

The spatially-resolved reflectance of turbid media has been studied experimentally and theoretically for 
many years. In particular, it has been shown that the measurement of the spatially resolved reflectance 
allows to determine the absorption and scattering properties of a turbid medium 1 " 9 . Such measurements 
have been applied for example in the medical field for oximetry 4 , photodynamic therapy 3,6,9 or glucose 
moni t oring 1 0 . Diffusion theuiy 1 " 7 , Muuie Cailu simulations 1 '^ 01 laiidom walk" have been used to 



describe theoretically the spatially-resolved reflectance. For source-detector separations larger than 
several transport mean free paths and high albedo, the diffusion equation gives accurate solutions for the 
reflectance, which depends in this case only on three parameters: the reduced scattering coefficient, the 
absorption coefficient and the index of refraction of the medium. 

The reflectance at shorter distances has been less systematically investigated, although it is of great 
interest. Indeed, measurements of the reflectance close to the source interrogates a smaller volume of the 
sample, compared to large distances measurements. Therefore, a local determination of the optical 
properties is possible. For example small volumes of biological tissues can be optically characterized 
during in vivo investigations 12,26 . 

At source-detector separations of approximately one transport mean free path, diffusion theory fails to 
describe accurately the spatially resolved reflectance. The problem of modeling the propagation of 
photons is then more complex than far from the source. Indeed, Bolt et al. 13 have clearly shown 
experimentally that the absorption and reduced scattering coefficients are not sufficient to describe 
accurately the reflectance close to the source. In particular, the phase function has to be taken into 
account. This has also been demonstrated by Kienle et al. 7,14 and Mourant et al 15 using Monte Carlo 
simulations. 

The first goal of this study is to assess the importance of the different moments of the phase function, in 



the spatially-resolved reflectance. In earlier studies 7 ' 13 , only the first moment of the phase function, also 
called anisotropy factor, has been considered to quantify the effect of the phase function close to the 
source. We will show that this analysis is not correct. Indeed, the first and second order moments of the 
phase function have to be taken into account, whereas moment of higher order can be neglected. 
The second goal of this study is to assess the validity of similarity relations, describing how the number 
of optical parameters can be reduced. We will show that, close t o the source, only three parameters 11 
beside the index of refraction are needed: the absorption coefficient the reduced scattering coefficient 
\i s ' and a third parameter y, defined as a function of the first two moments of the phase function. 
Therefore, only the additional parameter y is required compared to the case of large source-detector 
separations. 

The third goal is to describe how of these three parameters, p a , Ms* ^ Y can ^ determined from 
reflectance data. For this, the dependence the reflectance intensity, the slope of the log of the reflectance 
and the slope of the square of the reflectance is examined as a function of |i a , m s ' and y. Interesting 
properties are derived from this study and are be presented. 

The results of this theoretical study led us to the design of a probe measuring the scattering and 
absorption properties of tissues in vivo, with short source-detector separations 26,27 . Even if our main 
application target is the determination of the tissue optical properties, the results presented here have 
been obtained for a broad range of optical parameters, and can be applied to turbid media of various 
types. 

\ 



2. THEORY 
A* Geometry 

The geometry used in our simulations is the semi-infinite space. As described in Fig.l, we consider a 



light source and a detector placed normally to the surface and separated by a variable distance p. The 
spatially-resolved reflectance R(p) is defined by the power received by the detector per unit area and 
normalized by the source power. The numerical aperture of the source and the detector is 0.38 
(0 max =43.6°). This value has been chosen to match closely experiments where optical fibers are used to 
illuminate the medium and to collect the backscattered light 12 . 

B. Optical properties 

The medium has macroscopic ally homogeneous and isotropic optical properties: relative index of 
refraction n, absorption coefficient n a , scattering coefficient |i s and phase function p(9), where 8 is the 
scattering angle. From these coefficients we can define the total attenuation coefficient p t =^i s +ji a , the 
albedo a = p s /p t and the mean free path l/p t . 

The phase function can be expanded in a series of Legendre polynomials 16,22 ' 23 P n (G): 



(1) 



y g n are the rfi 1 order (Legendre) moment of the phase function: 



g n = 2rc£p n (0)p(0)sin0d0 



(2) 



The zero-order moment go is normalized to 1 for all phase functions. The first-order moment gj 
represents the mean cosine of the scattering angle 0. It is often called the anisotropy factor (or also 
asymmetry factor), and is usually simply noted g in optics. Thus, please note that g=gj . 
This parametrization of the phase function is interesting because, in most multiple scattering problems, 
only a limited number of moments need to be taken into account 16 . In particular we will examine the 
role of each moment on the reflectance, as a function of the source-detector separation. 



C. Choice of the phase function 
In this section we describe the different phase functions that we used for our study. Particularly, we 

define two phase functions, Pmhg anc * Pmpo the moments gj and g 2 of which, can be conveniently 

adjusted. They allow to approximate a broad range of real phase functions found in optics (from this 

point we most often drop the "0" dependence of the phase function for notation simplicity) 

Before introducing these phase function, we need m lerall tlieHeiiyuy-fi i t^ns i m i r jri i n.sK C om 
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which has been widely used for multiple-scattering problem. It allows to simulate highly forward 
scattering, and its mathematical form is simple: 

2 



m 



PHG( e )= 



1 



Who 



( 1+ gHG - 2 gHG COs9 ) 



(3) 



The moment of p HG are g n =gHG n where n is the order of the moment . The Henyey-Greenstein phase 
function is very convenient since its first moment can be easily adjusted. Nevertheless, fixing the value 
of gj determines also all the other moments. In order to investigate the effect of the second moment g 2 
independently of gj, we modified p H Q to a new phase function, we called Pmhg> hy adding a term 
proportionnal cos 2 9 (thus p^Q is a particular case of Pmhg)* 



3 2 

Pmhg( 0 )= ccPhg( 0 ) + ( 1 - (X )4^ cos 9 



a = 0..1 



(4) 



The weighting factor a guarantees the normalization of Pmhg* Th e term * n cos 9 gives a contribution to 
g 2 only, as for Rayleigh scattering 16 . Thus, Pmhg can he interpreted as an average between two types of 
scattering events: one with anisotropic scattering (e.g. originating from large particles compared to the 
wavelength) and one with quasi-isotropic scattering (originating from small particles compared to the 
wavelength). Phase functions with such a structure have been actually measured in biological 



tissues 
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The moments for Pmhg 2rc gi ve n as follows: 



Sl= «gHG>g2= a gHG 2 + 2 ( 1 - a )»g3 = «gHG 3 . g*= "gHC^ 



(5) 



Equ.(5) shows that the parameters g HG and a allow for the independent adjustment of gj and g 2 . 
Nevertheless the choice of gi and g 2 is limited by the condition that p(9)>0 for 0 < 9 < %. Due to this 
... A condition, > ,s.ome J . yalues .of gj^and^gg. ^hK^ ai;e^ n gt t 
possible with Pmhg, shown in Fig.2 In order to reach these values, we introduce a second phase 
function Pmpo constructed in a similar way to Pmhg : 



PMPc< e ) = ap pc (9)^l -a)^cos 2 8 a = °~ l 



(6) 



where ppc is a phase function made of power of cosines (corresponding to particular cases of the 
a w J expansion in Legendre polynomials of order N, see VandeHulst 10 ): 



1 N + 1 N 

PpcW = ^»^(l+cose) 



(7) 



Similarly to Pmhg> ^ e right term °f Pmpc proportioaal to cos 2 8 in Equ.(6) allows to adjust 
independently of gj . Indeed, the moments of Pmpc 316 tnen gi ve n by: 



gi= «gpc,>g2= agpc 2 + 2 ( 1 - <x )'g3= "gpc,. g4= «gpc 4 ' •• 



(8) 



The value of g PCn gpc depends on N: 



= N = (N-n+1) _ N(N- l)...(N-n+l) 

gpc > (N + 2) gpc n gpc »-'(N+n+l) (N + 2)...(N+n+l) 



(9) 



In Fig.2, the regions of possible g 2 are plotted for Pmhg PMPC ^ a function of g t . Values obtained 
from a broad set of Mie phase functions are plotted as well (relative refractive index from 0.9 to 2, and 
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size parameters from 1 to 25, from table 20 in Ref. 16). For high g] values, g 2 is generally characterized 
by g2^gi- Note that the region covered by the Henyey-Greenstein phase function ph G is only the line 
g 2 =gl 2 , which represents therefore very limited cases of Mie phase functions. 

Fig.2 clearly shows that Pmhg a U° ws to cover most of the g x and g 2 values from Mie theory (and also 
the particular case of Raleigh scattering characterized by gi=0, g2=0.1). Nevertheless, some cases are 
cannot be obtained by p M HG» but can bc r eached using Pmpc- Generally, -the cases covered by Pmhg 



also covered by Pmpc- For such cases, the use of Pmhg m ^ PMPC 1S complementary. Indeed, for identical 
g! and g 2 values, the higher oreder moments g n (n>2) given by Pmpc ' m significantly lower than the 
ones given by Pmhg* 1S shown in Fig.3, for the cases gx=0.9, g2=0-Bl and gj=0.5, g 2 =0.25. 
Therefore, the use of either Pmhg or Pmpc is motivated depending on the choice of the moments of high 
order (> 2). In particular, we will use Pmhg and PMPC to assess the role of moment of order higher than 
two on the reflectance. 



D* Similarity relations and reduced coefficients 

Consider a turbid medium with optical properties (p a , p s , p(0)). Is it possible to find a different set of 
optical properties (|i a *, p s *, p(9)*) that gives the same reflectance profile? If the answer is positive, the 
relations between these two sets of parameters are called the similarity relations. This concept, 
introduced by Van de Hulst 16 , has been well discussed by Wyman et al. 22 * 23 The similarity relations are 
interesting because they allow for the reduction of the number of parameters necessary to describe the 
reflectance. 

Wyman et al 23 have shown that if the radiance L(r, s) can be expressed as a finite spherical harmonics 
Y™(0, (p) expansion (cp is the azimuthal angle): 

N n 

L(r,s) =XI a mn (r)Y>,<p) (10) 

n - 0m = -n 
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then it is possible to define N+l similarity relations: 



O 



in 

5 == 



Ma* = Ma- 



H s *(l-g„*) = M s d-g„) 



n= 1.....N 



(11) 
(12) 



where g n are the moments of order n of the phase function, as defined by Equ.(2). 



If N=«>, the exact phase function has to be used to obtain the correct radiance distribution. 

N=l is the well known case of the diffusion approximation, where the radiance is only linearly 

anisotropic. Only two similarity relations are necessary (called first order similarity relations): 



Ma* = M a 



Ms*a-gi*) = M s a-gi) 



(13) 



(14) 



In this case, it is useful to define the following reduced coefficients: 

- the reduced scattering coefficient: m'=H s (l-gi), 

- the reduced total attenuation coefficient: \x^~\x t (l-gi), 

- the reduced albedo a*=|i s 7(i t \ 

- the transport or reduced mean free path mfp'~l/\i t y 

When the diffusion approximation holds, media with identical reduced coefficients lead to identical 
radiance distribution. For this reason, n s ' and are commonly used to characterize optically thick 
turbid media. 

The diffusion approximation (N=l) is generally valid for source detector separation larger than several 
transport mean free paths. For shorter distances, the case N=2 (radiance at most quadratically 
anisotropic) should be considered, then N=3 and so on. 
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For the case N=2, a third similarity relation should be added to Equ.(13) and Equ.(14): 



Ms*(l-g2*) = M s (l-g 2 ) 



(15) 



Combining Equ.(14) and Equ.(15), the third similarity relations can be rewritten as: 



l-g 2 l-g 2 * 



=5? 



(16) 



Therefore in case N=2 (radiance at most quadratically anisotropic), three parameters have to be 
considered: p s \ and y = (l-g^O-gi). The validity of these similarity relations will be assessed by 
simulations. 

The y value depends on the balance between the first and second moment of the phase function. 
Examples of y values, computed from Mie scattering, are given in Fig.4 for monodisperse spheres. In 
this figure, y is plotted versus the parameter diameter/wavelength, for three different index of refraction 
ratios n = n sphere /n medium : 1.05, 1.1 and 1.2. For very small particles, y decreases if the size of the 
particle decrease and has a limit value of 0.9 for Rayleigh scattering. In such a case, the first moment 
gl=0, but the second moment is g2=0.1. For poydisperse scatterers, the parameter y is also sensitive to 
size distribution. In the particular case of two populations of scatterers, one of very small size (Rayleigh 
scattering) and one with large particles, the average phase function can be approximated by Pmhg or 
Pmpc- From Equ.(5) and Equ.(8), it can be easily shown that y is very essentially sensitive to the ratio of 
the concentration between large particles and small particles (this ratio is equal to cc/(l— a)). 



E. Monte Carlo model 

The Monte Carlo method is well known to give accurate solutions for the light propagation in turbid 
media, if interference effects are negligible. The Monte Carlo code we used has been extensively tested 



and compared with other codes. 
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The source and detector characteristics follow the description given in Fig.l. Boundary conditions are 
taken into account using Fresnel and Snell laws for each photon reaching the surface. 
Taking advantage of the semi-infinite geometry and cylindrical symmetry, scaling procedures can be 
used to limit the number of simulations. The idea is simply to express the source-detector separation by 
the unitless quantity p - |a s ' . Indeed, for a given reduced albedo a\ the unitless quantities 



-i 



(|i s ') "R(p|X s ', a') computed trom a single simulation (for exampl e chuw iii g - 1 mm" , d - 0 . 99) , 

allows to derive the reflectance R(p) for any values (but constant reduced albedo a'). This is 

explained as follows. Consider the spatially resolved reflectance R(p, p s \ a') obtained with a medium of 

coefficient n s ' and reduced albedo a'. The reflectance R(p, k ji s \ a')> obtained with a medium of the 

same albedo but with a scaled coefficient k p s \ is simply related to the previous one by rescaling p by kp 

28 

and by renormalizing the reflectance (because of the rescaling of the area of the detectors): 



R(p, kji,', a") = -^R(kp, n s \ a') 

V 



(17) 



—2 

We present our results of the reflectance using such unitless quantities ((i s *) R(p • \i s \ a') or 
p 2 R(p • \i\ a' ) , in order to obtain the most general results. Similarly, the dependence of slope of the log 
of the reflectance as a function of n s ' or p can be derived the unitless quantity p^lnR(p • |A S \ a') . 



3. SIMULATION RESULTS 

A. Effect of the moments g n on the reflectance 
First we examine the effect of the different moments g n on the reflectance. Results for the reduced 
albedo a' =0.99 are plotted here. Such albedo are found in biological tissue for red and near-infrared 
wavelengths. The effect of the absorption on the reflectance is examined later (section C and section D). 
A relative refractive index n=l .4 of the medium is considered. When indicated, the results obtained with 
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n=1.4 are compared to the case n=l. 

The effect of gj is first examined. For this, g 2 is held constant and gj is varied. In Fig.5, the reflectance is 
computed with pmhg with gj=0.9, g 2 =0.81 and p MHG w i th gi=0.81, g 2 =0.81. Fig.5 shows clearly that 
the influence of g x is important at distances around pp s '«l. Note that such differences are much larger 
than if the Henyey-Greenstein phase function was used for this test, as it has been performed by some 
authors This is due to the fact that when p HG is used, g 2 is not held constant when g } is changed, but 
varies as g x 2 . In such a case, as it will be clarified by the similarity relations shown in section B, the 
effect of gi is underestimated. 

The role of the second moment g 2 is examined in Fig. 6.a and Fig.6.b (Fig.6.b is a close view of Fig. 6.a, 
for 0<pp s '<2). All the reflectance curves are obtained with phase functions of identical gj = 0.9 but 
different g 2 values: 0.75, 0.8 1, 0.9. p MH c was used for &2 = 0 9 ^ d °- 81 (Phg=Pmhg in this case )> mA - 
PMPC in the case of S2 = °- 75 * As shown in Fig.2, these g 2 values correspond, to extreme cases found 
when considering Mie phase functions with high g v The solution of the diffusion equation given by 
J? Kienle et al. 30 is also presented in these figures, in order to quantify the limit of validity of the diffusion 
theory. The angular distribution of the reflectance was assumed as Lambertian 7 . 

Fig. 6.a clearly shows that all the reflectance curves converge to the diffusion solution for pp s * > 10 . It 
demonstrates that for such distances the diffusion approximation is valid. Nevertheless the divergence 
between the diffusion solution and simulations depends highly on the second moment g 2 . For example, 
in the case of g 2 = 0.9, the reflectance is approximately 25% smaller than the diffusion solution at 
PH S '=2. In contrast, the diffusion solution matches closely the reflectance obtained with p HG , where 
differences smaller than 5% are found for distances as close as pn s ' = 0.5. The importance of g 2 is 
clearly demonstrated in Fig.6.b g 2 can induce differences in the reflectance up to 30% at distances 
0-5<PM s ' < 2.The role of g 2 is found here to be similar to that of g^ shown previously. Therefore, the 
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combination of both gi and g 2 is responsible for the reflectance curve shape at pu s '~l . 
One could suspect that the differences found in Fig.5 or Fig.6.b are not only due to differences in gj or 
g 2 , but also to differences between the moments of higher order. To study the role of the moments g n 
(n>2), two sets of phase functions with identical g x and g 2 , but different higher moments are used. Two 
cases were investigated: g^O.5, g 2 = 0.25 and gi=0.9 and g 2 =0.81. For this test, Pmhg PMPC 
used, because their moments g n exhibit important difference s lor il2:3 (see Fig.39j. Such difference s are 



close to maximum values that can be found between different Mie phase functions. Fig. 7.a and Fig.7.b 
show that the reflectances curves computed using these two phase functions are very close for pu s '>0.5 
(differences smaller than 5% in Fig. 7.a and 15% in Fig.7.b). For smaller distances, the differences 
become important (>50% for pu s '<0.1). Compared to the effect of g x in Fig.5 or g 2 found in Fig.6.b, the 
effect of g n with n>2 is weak for the reflectance at distances pus'>0.5, whereas the g n values (n>2) have 
definitely to be taken into account to compute the correct reflectance at smaller distances. 



B. Similarity relations 

The fact that only the two first moments play a significant role for the reflectance at distances larger than 
0.5 mfp' suggests that the second order similarity relation may be valid. In such a case, the similarity 
relations given by Equ.(13), Equ.(14) and Equ.(16) indicate that the reflectance should depend only on 
the parameters: Y=(l-g 2 V(l-giX H s ' and Ma- To demonstrate the validity of these similarity relations we 
report in Fig. 8.a the reflectance computed with three phase functions characterized by identical 
y= 125, but with gj = 0.2, 0,5, 0.9 and g 2 = 0, 0.375, 0.875 respectively. 

Fig. 8.a shows that these similarity relations are satisfied within only 2% error margin for pu s *> 0.5, in 
the case n = 1.0. For the case n=1.4, slightly higher differences are found: <10% for p|i s *> 0.5. Note that 
these differences are much lower that differences found in Fig.5 and Fig.6.b, where g x and g 2 where 
varied independently. 



14 



*=5 f 



Using these similarity relations, isotropic scattering (g n =0 for all n, 7 = 1.0) can also be approximated by 
a phase function characterized by 7= 1.0, even with high gj. This is illustrated in Fig.8.b, where the 
reflectance produced by isotropic scattering (g n =0) f and Pmhg gi = g2=0-9 are compared. As before, 
the similarity relations are found efficient. The differences found are of the same order than in Fig. 8. a 
In view to these results, the hypothesis of quadratically anisotropic radiance for [i s '> 0.5, on which is 
based the simila r ity relations wc used, appear to be well satisfied, hor n = 1. 4 , the mismatch or retractive 
index modified the angular distribution of the radiance, which explains that the similarity relations are 
less perfecdy satisfied for n=1.4 compared to the case of n=1.0 . 

In summary, Fig. 8.a and Fig.8.b clearly show that the effect of the phase function can be quantified by 
only the parameter 7 only, for distances p|i s *>0.5, constant fi s ' and a*. For shorter distances, the 
moments of order higher than two must be taken into account, which would involve other similarity 
relations to be considered. 

C. Effect of the albedo on the reflectance. 

We study in this section the effect of the albedo on the reflectance at short distance. Increasing the 

absorption induces changes in the intensity and in the slope of the reflectance. To examine this effect 
quantitatively, we derived from the reflectance the two unitless parameters p 2 R(pn s \ a') and 
p9plnR(p|i s \ a*). For p=l mm, these parameters represent respectively the intensity of the reflectance 
and the slope of the log e (noted In) of the reflectance, measured at 1 mm. Scaling relationships allow to 
derive easily these latter quantities for other distances. 

The parameters p 2 R(p|n s \ a') and lp3 p lnR(p|i s \ a')l are plotted in Fig.9 for reduced albedo from a'=l to 
a'=0.83 and for distance pps' from 0.28 to 5, and for two different phase function with the same first 
moment gj = 0.916 (p HG with g^.916, 7=1.92 and p^ with g^O.916, 7=2.23). 

For distance \i s y p larger than approximately 5, the difference induce by these phase function are 
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negligible. This is in agreement with the results presented in section A. At shorter distance significant 
differences can be found, depending on the phase function used. As shown in the previous section, y is 
the important parameter of the phase function to consider, if the distance is larger than approximately 
Ms 'p = 0.5. 

Fig.9 shows that, for a fixed phase function and a fixed distance p, the two parameters p R(pn s ') and 



p8pliiR(pH s T ) of the ate related uniquely to p[i s J and a' . Therefore, tor a given value 01 p, p~K(pn s 'j an<F 

p3plnR(pH s ') allow to determine the parameters a' and \x s " (or equivalently |i a and ji s ') of the medium 

with known phase function. If the phase function is not known, the determination of a' and \i s y require 

the measurements of p^Cpp/), p9prnR(p|a s *) at a given distance p and at least an additional parameter. 

This latter one could be for example p 2 R(pn s ') or pd p lnR(pn s ') measured at a second distance p. This 

point is discussed further in the section 4. 
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D. Approximate form of the reflectance close to the source. 

We present in this section empirical properties of the reflectance close to the source, obtained by 

simulations. They may be used to determine the optical properties form reflectance data. For reduced 
albedo a'>0\9 (corresponding to tissue reduced albedo found for red and near-infrared wavelengths) and 
pji s '<4, we found empirically that the reflectance can be expressed in the form: 



R(p, n a , jl \ y) = [A(p, \i s \ y) + B(|a s \ H a )] 



(18) 



where A and B are functions that are described in Fig. lO.a and Fig.lO.b. The function A depends on the 
scattering properties (i.e. p s ' and y) but not on the absorption. In contrast, the function B depends on the 
absorption but not on the phase function (i.e. y)- Equ.(18) is interesting because it shows that it is 
possible to uncouple the effect of the phase function and of the absorption. We illustrate this 
approximation in two steps, by deriving two property of Equ.(18). 
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First, the slope of the square root of the reflectance - — Jr does not depend on the absorption coefficient 

dp 



(19) 
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The dimensionless quantity p"^-(jR) is plotted in Fig. lO.a as a function of pji s ' for y =1 , 1 .9 and 2.5 

dp 

and icduced albedo a'-l , 0.95 and 0.9. - 
Fig. lO.a shows that p 2 ^-{jR) depends highly on p|i s ' and y In contrast, the dependence on a' is 



dp 



almost negligible. If y is known, jj s * can be simply monitored by the parameter p ^(a/R), 
independently of the variation of p a . Moreover, the inversion procedure, i.e determining p s * from the 
H reflectance curve is simple using p ^— (JR.) , because this parameter can be well approximated by 

u p 

H polynomial functions (of order 4 or 5) (for the range 0.4<pp s '< 4). 

W 

Ul The second properties derived from Equ.(18) is as follows: 



yj 



JR(P,^, \i s \y) - JR(p, ^ a = o, p/,y) = #B(^ a , ^i;)-B(p a = o, ^) 



(20) 



3 The difference jR(p, Ji a , H S \Y) - ^(p, |i a = 0, ^ s ', y) can be transformed in dimensionless quantity 
B by multiplying it by p. This dimensionless difference illustrated in Fig.lO.b for the same three phase 
functions that we used in Fig. lO.a and for a'=l to 0.83. Fig.lO.b shows that the influence of the phase 



function is weak in the quantity p(jR(p\i s \ a') - jR(p\x s \ a'= 1)). Moreover the relation found in 
Fig.lO.b is very close to a linear function. Therefore, the further approximation can be made: 



(VR(pm\ a') - VR(pn s \a*= 0) = ^'f (a 1 ) 



(21) 



where f(a') is a function depending only on a'. For known p s \ the function f(a') allows the 
determination of a relative absorption change Ap a =p a - p^, from a known value p ao . Fig.lO.b illustrates 
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the case p ao =0, but any other vaJue of |J ao is possible. The interesting point is that f(a') does not depend 
on the phase function. 

The two properties given by Equ.(19) and Equ.(20), proven by simulations examples in Fig. lO.a and 
Fig.lO.b, respectively, confirm the form of the reflectance (close to the source) given by Equ.(18). 



4. DISCUSSIONS AND CONCLUSIONS 

The study of the spatially resolved reflectance is easier if distinct ranges of distances are considered, 
corresponding to different regimes of the photon propagation. Each regime is related to a specific 
m angular distribution of the radiance, which varies from a highly anisotropic distribution (close to the 
rf source) to a quasi-isotropic distribution (far from the source) 31 . In this discussion we consider three 
[j } different regimes . 

'j\ The first regime occurs where the diffusion approximation holds and corresponds to a linearly 
= anisotropic radiance. Its range of validity can be determined indirectly by testing the validity of the first 
^ order similarity relations (^i a =cte, \x s y = cte). Simulations performed with a set of different phase 
Ji functions show that their effect on the reflectance is less than 5% for distances defined by p|i s '>10 as 
long as ^i s * and (J a are kept constant (Fig. 6.a). The different phase functions we used for this test 
approximate phase functions that can be obtained from Mie theory or phase functions that have been 
measured in biological tissues. 

This limit of the diffusion regime is comparable to other estimations obtained when studying other 
quantities such as the spatially resolved transmittance 24 or the radiance angular distribution inside a 
turbid medium 32 . Nevertheless, this estimation is quite conservative. As demonstrated by other 
authors 5,30,33 , solutions of the diffusion equation can still be accurate much closer to the source. In 
particular, in agreement with Kienle et al 30 , we found that the diffusion solution is very close to 
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simulations down to distances PMs'-l, if Henyey-Greenstein phase functions with high anisotropy 
i factors are used. However, this not true for other phase functions. Thus, diffusion theory should be used 
with care for distance pju^lO, if the actual phase function is not known. 

The second regime is delimited approximately by 0.5<pp s '<5 range, and corresponds to the region 
where the radiance is approximately quadratically anisotropic. In this case, a third similarity relation can 
be used: the reflectance depends on ttic parameter 7 ° U §2 ) / ( 1 £ 1 )> which must be ke P l constant in 
addition to n s ' and p a . Therefore, both the first moment and the second moment g 2 of the phase 
function must be taken into account. In this regime, the effect of moments of higher order is weak on the 
reflectance, compared to the effect of g] and g 2 . 
Jl! The third region corresponds to distances p|i s '<0.5. The anisotropy of the radiance becomes high, and 
p more moments must be taken into account. In such a case, computations with accurate phase functions 
Ul should be performed. 

^ Our analysis also demonstrated that it is not possible to determine the value gj alone by the 
^ measurement of reflectance close to the source, if no other parameter of the phase function is known. In 
rj view of this, we believe that the errors found by Jones and Yamada 29 , in the measurement g t by such a 
m technique, are explained by the fact that they neglected the effect of the second and higher moments in 
their theoretical analysis. It also shows that testing the effect of g x using the Henyey-Greenstein phase 
function alone is not correct. Indeed, in this case, both g! and g 2 =gi 2 are varied simultaneously. For g\ 
between 0.8 and 1, y varies then only from 1.8 to 2.0. This explains that the influence of gi on the 
reflectance has been reported 7 * 34 to be weak for g^O.8. However, as shown in Fig.5, the effect of gj can 
be much larger if g 2 remains constant. 

The three coefficients jj s \ ji a , and y can be determined if measurements are performed in the diffusive 
regime and close to the source (pn s '>0.5 rrifp') on homogenous samples. [i s \ ji a can be obtained by 
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fitting the diffusion solution to the measurement at large distance 5,7 * 33 ' in order to avoid the dependence 
of y. The y value can be then determined by comparison between the reflectance measurement close to 
the source (pp s ~l) and a set of simulations performed with various y values. 

If the measurement is performed only close to the source, the determination of the p s \ p a , and y are 
possible by fitting simulations to the measured curve, but with a lower degree of precision. Note that 
non - linear regression could also be performed using the method proposed rec ently by Kicnlcct al.-^ F 



which consists in using ad hoc analytical functions fitting the reflectance simulations and scaling 
relationships. 

At short distance, the problem of computing the optical properties is facilitated if y or p s ' is known. If y 
is known, ji s * and |i a can be derived from the reflectance intensity R(p) and the slope of the logarithm of 
the reflectance 3 p InR(p) at one distance p close to the source. From Fig.9, the precision of such 
determination can be estimated. For an optical distance of 1 mfp\ a precision better than ±1.5% is 
achievable on the determination of |i s ' if these two parameters are measured with an accuracy of ±1%. 
The accuracy on p s ' is weakly dependent on p a and p s \ for typical tissues properties (0.5< p s *<2 mm" 1 , 
[i a <0.2 mm" 1 ). In contrast, the accuracy achievable on p a is highly dependent of the |a a value. At a 
distance of 1 mfp' and p a =0.1 mm -1 , the error on |i a is approximately ±5%, taking into account 
uncertainties of ±1% on 9 p lnR(p) and R(p). For p a around 0.02 mm" 1 , the errors increases to 
approximately ±20%. 

If p s * is known, but not the phase function, relative measurements of |i a can be still achieved using the 
empirical property demonstrated in section D (Equ.(21)). This feature is very interesting for absorption 
monitoring in tissue for example. 

In conclusion we showed that the spatially-resolved reflectance region at distances close to one transport 
mean free path can be accurately predicted with only three parameters: p a , p s \ 7= (l-g2V(l-gi)* These 
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results were applied for the determination of [i a and |i s * of biological tissues 26 , using spatially-resolved 
measurements at distances 0.5<p<1.5 mm, corresponding in tissue typically to distances of 
approximately one transport mean free path. 
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FIGURE CAPTIONS 

Fig.l Geometry of the study. Photons are emitted at p=0 isotropically in a solid angle defined by the 
^gte e max- They are detected at variable distances p in the same solid angle defined by 9 max . 

Fig. 2 Possible values of g 2 as a function of g, for Mie phase functions and Pmhg m ^ Pmpc- The 






region covered bv the our is renresented bv the line gorp,^ The line renrpcent the nh»^ 






function formed by isotropic scattering added of a purely forward peak 16 . The dots represent a 






broad choice of Mie phase functions (relative refractive index from 0.9 to 2, and size 






parameters from 1 to 25, from table 20 in Ref. 16) 




Fig.3 


Values of the moment g n for Pmpc and p MH G- Th c values of the moment g G , gj and g 2 are 






identical for both phase functions. Two cases are shown: g l =0.9, g 2 =0.81 (upper curves) and g l 






=0.5, g 2 =0.25 (lower curves). 


o 
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diameter/wavelength, for different index of refraction ratios n = n sphere /n medium . 




Fig.5 


Effect of the first moment g=g] for constant second moment g 2 . 




Fig. 6.a 


Role of the second moment on the reflectance, for a fixed first moment g } = 0.9. 




Fig.6.b 


Role of the second moment on the reflectance, for a fixed first moment gj = 0.9. This graph is a 






close view of Fig. 6.a 




Fig. 7.a 


Effect of the third and higher order moments on the reflectance. The first two moments are 






fixed: g 1= 0.5, g 2 =0.25. 




Fig.7.b 


Effect of the third and higher order moments on the reflectance. The first two moments are 
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fixed: gj= 0.9, g 2 =0.81 

Fig. 8. a Illustration of the second order similarity relations. All phase function have identical parameter 
T=(l- g2 )/(l- gl ) = L25 

Fig.8.b Illustration of the second order similarity relations. The isotropic scattering can approximated 

by a phase function with high gl if gi=g2- In such a case y = (l-g2)/(l-gi) = 1- 
Hg.y Relationship between the parameters p 2 k(pu ' ) and pdJkikipu' ) and the reduced scattering 



coefficient jj s ' and the reduced albedo a*, for p HG with g=0.9 (7^=1.9). Mismatched boundary 
condition n=1.4. 

^ Fig. lO.aPlot of the parameter p 2 -^~(jR) for different phase functions and different reduced albedo. The 

y ■■ dp 

m= following phase functions were used: Pmhg w ^ Si = 0-9, g2 = 0.9 Cpl.O), Phg w ^ gi=0-9,- 

W g 2 =0.81 (7^=1.8) and Pmpc with gj=0.9, g 2 =0.75 (y=2.5). Mismatched refractive index n=1.4. 

HJ Fig.lO.b Plot of p( A /R(pji s , } a') - jR(p\i^ 9 a'= 1)) for different phase function and reduced albedo a'. 
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ABSTRACT 

Local and superficial near infrared (NIR) optical property characterization of turbid biological tissues can be 
achieved by measuring spatially resolved diffuse reflectance at small source-detector separations (<1.4 mm). 
However, under these conditions, the inverse problem, i.e. calculating localized absorption and the reduced 
scattering coefficients, is necessarily sensitive to the scattering phase function. This effect can be minimized 
if a new parameter of the phase function y which depends on the firs t and second moments ot the phas^= 
function, is known. If y is unknown, an estimation of this parameter can be obtained by the measurement, 
but the uncertainty on the absorption coefficient is increased. A spatially-resolved reflectance probe 
employing multiple detector fibers (0.3 to 1.4 mm from the source) is described. Monte Carlo simulations are 
*p used to determine y the reduced scattering and absorption coefficients from reflectance data. Probe 
P performance is assessed by measurements on optical phantoms, the optical properties of which was 
R measured by other techniques (frequency domain photon migration (FDPM) and spatially-resolved 
Hj transmittance) and reasonable agreements are found. However, the spatially resolved probes shows optimum 
measurement sensitivity in the measurement volume immediately beneath the probe, while FDPM typically 
N= samples larger regions of tissues. In vivo measurements performed intraoperatively on human skull and brain 
are reported for 4 NIR wavelength (674, 811, 849, 956 nm), using both methods. Optical property values 
J for human skull, white matter, scar tissue, optic nerve and tumors are reported, which show distinct 
" absorption and scattering differences between structures and a clear dependence on the phase function 
parameter y. 

Keywords: turbid media, tissue optics, optical biopsy, neural tissue, absorption coefficient, scattering 
coefficient. 
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1. INTRODUCTION 

' Probing the optical properties of biological tissues has a major impact in several medical applications for 
diagnosis and therapy. For example, the knowledge of these properties is necessary for optimizing 
techniques such as near infrared spectroscopy or photodynamic therapy. The scattering and absorption 
characteristics of many different kinds of tissues have been reported in the literature 1 . However they have 
been mostly measured in vitro. Because of unavoidable alterations of excised samples, such as blood 
drainage, structural alterations, and temperature changes, these values are questionable and in vivo 
measurements are preferable. 

The measurement of optical properties performed in vivo could be also used as a diagnostic tool, which is 
complementary to other optical biopsy techniques, e.g. tissue autofluorescence. Indeed, light scattering and 
Cn absorption can provide information both on tissue structure and chromophore content, features which can be 
Mused to distinguish between normal tissues, malignant lesions and other pathologies. For example 
Ld hemoglobin and water content have been found to be significantly different in normal and cancerous 
y=i tissues 23 . Differentiation between normal and malignant bladder tissues 4 have found to be possible from the 
J* elastic scattering and absorption properties. 

El Different methods have already been proposed to quantitatively determine the absorption and reduced 
H scattering coefficients in vivo, using spatially- 5 " 7 and/or temporally-resolved measurements 2,3 * 8,9 . Besides 
organs such as breast or brain which can be transilluminated, measurements of thick tissues can be made in 
reflectance geometry. The case where the source-detector separations are larger than several transport mean 
free paths, corresponding typically to distances larger than 5 mm for biological tissues, has been extensively 
studied theoretically and experimentally. Diffusion theory or Monte Carlo simulations have been commonly 
used to relate the measured light intensity to the optical coefficients. In particular, optimized source-detector 
separations have been calculated by different authors 6,10 in order to ensure the best sensitivity to absorption 
and scattering coefficients from spatially-resolved reflectance data. For typical turbid tissues, optimal 
determination of both absorption and scattering properties requires reflectance measurement at small and large 
distances. For example, measurements were performed by Patterson and co-workers using distances from 1 
to 10 mm, and by Bays and co-workers using distances from 3.6 to 15 mm or 2 to 14 mm. 



Nevertheless, all these studies consider the turbid medium as homogeneous, or made of homogeneous 
» layers. Our approach is different. We wish to optically characterize a small volume of tissue, on the order of 
a few cubic millimeters, possibly distinct from the surrounding tissues. Therefore, our goal is to differentiate 
a small tissue heterogeneity instead of determining average optical properties of a large volume of tissue. To 
achieve this we chose to perform spatially resolved reflectance measurements with only small source-detector 
separations, from 0.3 to 1.4 mm, even if these distances are not optimal for absorption determination. 
Mourant et al 4 U have shown that the absorption coefficient can be estimated with a measurement at a single 
distance of approximately 1.7 mm, assuming the scattering coefficient to be in a certain range. They have 
also shown that measurements at a single shorter distance (approximately 0.3- 0.4 mm) allows monitoring of 
spectroscopic variations in the scattering properties of tissues. In this work, we address a general case where 

uiboth tissue scattering and absorption properties are simultaneously estimated, from measurements at 

p distances smaller than 1.4 mm. 

M A previous study we performed using Monte Carlo simulations gives the theoretical framework for the 
interpretation of the measured profile in terms of absorption and scattering properties 1213 . The role of the 
s4J phase function at short source detector separations were carefully studied. In particular, we showed that, 
t! beside the absorption and reduced scattering coefficients, a parameter depending on the first and second 

IssP 

s moment of the phase function must be taken into account, for source detector separation ranging from 0.3 to 

SI 

*fi 5 transport mean free path. This analysis is more complete compared to previous work which considered 



only the first moment of the phase function (anisotropy factor). 

We first describe the probe design and the Monte Carlo model developed to simulate the measured profile. 
Second the average tissue volume probed, and sensitivity to medium boundaries are discussed from 
experiments and simulations. Third we describe how the scattering and absorption properties of tissues can 
be deduced from reflectance data. The accuracy of the proposed procedure is demonstrated on phantoms, the 
optical properties of which are measured by others techniques. Finally measurements of human brain and 
cervix tissues obtained in vivo are presented and discussed. These clinical measurements were performed in 
parallel with a complementary method, frequency domain photon migration (FDPM) 2 * 3 , which probes a 
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larger tissue volume. Optical properties from FDPM measurements are compared with values obtained by the 
spatially-resolved method described here, using small source-detector separations. 



2. MATERIALS AND METHODS 



A. Definitions 



The spatially resolved reflectance is denoted R(p) where p is the source-detector separation. It is defined by 
the power received by a detector per unit area for a source of power unity. In our measurements p ranges 
between 0.3 and 1.4 mm. 

The optical properties of tissues are the average index of refraction n of the medium, the absorption 
^coefficient p. a , the scattering coefficient n s? and the phase function p(0) where 9 is the scattering angle. The 
El phase function is the density probability function for 0. We consider the index of refraction of tissues 14 as a 

constant n=1.4. 

Hht is also useful to define the reduced scattering coefficient |i s *= jx/l-g) and the transport mean free path mfp* 

4J 

£ = (iV+mV 1 where g is called the anisotropy factor and is defined as the average of cos0. Generally the 
□ reduced scattering coefficient |X $ ' and the absorption coefficient |i a are used to characterize optically thick 
SI tissue. Indeed, the light fluence rate depends only on |i s ' and (i a at distances of several transport mean free 
Cri paths (typically p>5mm) from the source (diffusion approximation). Therefore the use of |i/ and \i a is a 
natural choice if measurements are performed at such distances. As we want to make measurements at closer 
distances, in the range of one transport mean free path, we expect that some parameters of the phase function 
have to be taken into account. This theoretical problem was fully studied with Monte Carlo simulations and 
reported in reference 12. We consider below only the main implications of this work. 

For distances between than 0.3-5 [mfp'], we found that the reflectance curve depends on ji a , |X $ ' and a third 
parameter y = (l-g 2 V(l-gi)» where g l and g 2 are respectively the first and second moment of the phase 
function. The parameter y is derived from the second order similarity relations derived by Wyman et al 15 , 
which are valid for a quadratically anisotropy radiance. For comparison, note that the diffusion 
approximation and first order similarity relation correspond to a linearly anisotropy radiance. 



Generally, the n* moment g n is defined as 16 : 



g n -27tfp„(0)p(e)sinede 

0 



(1) 



where P n is the Legendre polynomial of order n. Note that the first moment g x is the anisotropy factor g. The 
role of the parameter Y=(l-g 2 )/(l-g t ) implies that the anisotropy factor g (=g,) alone is not sufficient to predict 
correctly the reflectance curve close to the source. 

From this analysis, each tissue can be potentially characterized by three parameters: |l a , and y. However, 
due to the restricted range of the source-detector distances we want to use, the simultaneous determination of 
H a , m' and y is not always possible with high degree of accuracy. The achievable accuracy depends on the 
optical properties themselves and on the experimental uncertainties. This problem is addressed in the results 
: section 3. for both phantom and tissue measurements. 
R The parameter y may give interesting information about the tissue structure. Indeed, as measured by several 
H3 authors 17 " 20 , tissue phase function can be seen as a sum of a highly forward phase function p HA ( 6 )> due to 
^ large particles, plus a low anisotropy phase function p^B), due to small particles. 



P^ucO) =(l-a) PHA (9) + a Pla (6) 



(2) 



The coefficient a is introduced to guaranteed the normalization of p ussuc (9). 

The first term p^G) has been often fitted to the Henyey-Greenstein phase function: 



Phg(9) = 



1-gHo 



(l + gHG- 2 gHG COS G) 



(3) 



The moments of the Henyey-Greenstein phase function are given by g„ = g HG n (n>0). 

The second term Pla(9) can be interpreted by Raleigh scattering, valid for very small scatterers compared to 
the wavelength. It has been approached by a purely isotropic term 17 or an Henyey-Greenstein phase function 
with a low negative g value 18 . We propose to use the exact Rayleigh phase function 16 : 

6 



p Raytegll (e)=^(i+cos 2 e) (4) 

The moments of the Rayleigh phase function are: gj= 0, g 2 = 0.1, g 3 = 0, g 4 = 0, ... 

Rayleigh scattering give no contribution to first moment g lf but only to the second moment g 2 . The moments 
of the tissue phase function p dssue (6) given by Equ.(2) are therefore: 

g, = (i-cQg™, g, = (i-a)g H « a + Q-i« (S) 

Therefore, we see that the parameter y is influenced by the relative concentration of Rayleigh scatterer a, 
which should depends on the tissue structure. 

Published phase function data suggest possible values of y for biological tissues 17 20 . They have been 
— performed by goniometric experiments of thin samples. Many artifacts can affect these measurements, such 
J^as the tissue preparation and the tissue thickness, and these results should be used with caution. The phase 
p; function reported by Jacques et al. 17 for human dermis at 633 nm leads to y = 1 .4. The phase functions of 
U] white and grey matter at k =750 nm measured by van der Zee et al. 19 give values of approximately 1.5. 

Therefore we performed p MHG simulations with y = 1.5 as a starting point. As will be discussed, fitting our 
M experimental data to simulations performed with different values of y estimation, permits reasonable 
q parameter estimation. 

2? B . Experimental setup 

The probe used for the measurement of the spatially resolved reflectance is described in Fig.l. It is a linear 
array of optical fibers (core diameter of 200 Jim, N.A. = 0.37 in air). Two source fibers can be used to 
illuminate the tissue. They are disposed symmetrically in regard to the collecting fibers. If the sample is 
homogeneous, the reflectance curve is identical with either illuminating fiber. Therefore, comparing the two 
curves tests the heterogeneity of the investigated tissue region or detects obstructions, beneath the 
iUuminating fibers. If the two curves are close (typically differences less than 10%), the measurement is 
validated and the average of the two curves is calculated. 

The illuminating fibers are slid in small stainless steel tubes in order to avoid direct light coupling with the 
collecting fibers. The coupling between each collecting fiber has been experimentally measured and found to 



be less than 2%. The fiber array is set in a stainless steel tube of 2.2 mm diameter and 20 cm long. The tube 
■ is filled with an optically clear adhesive. The probe is rigid, which allows for easier handling by the 
physician, during surgery for example. The whole probe can be sterilized. 

The experimental setup is shown in Fig.2. An optical switch (Dicon, model GP700) is used to select the 
illuminating fiber from different sources. For the brain measurements, four laser diodes emitting at 674 nm, 
81 1 nm, 849 nm and 956 nm were used (SDL, Inc. models 7421, 5420, 5421 and 6321, respectively). Two 
other laser diodes were used for the phantom measurements, emitting at b /3 nm ana at il& nm (ILEE LDA 
2011 and 1805 respectively). The six fibers used to collect the backscattered light are imaged on a linear 
Charge-Coupled-Device (CCD) (Hamamatsu S3921). The signal is digitized by a 12 bit A/D card. Only one 
measurement, which takes approximately 0.1 s, is then needed to measure simultaneously the intensity 
m collected by the six fibers. The entire system is controlled by a personal computer. 

2 Transmission differences between each fiber are corrected using a measurement on a turbid phantom 
bj illuminated uniformly. Immediately after each reflectance measurement, a measurement of the background 
m light is automatically performed and then subtracted from the reflectance signal. To minimized the 
f background light, a long pass filter (X > 650 nm) is put between the end of the bundle and the CCD. Even 
K during open surgery where the ambient light is substantial, the measured background was less than 5% of the 
==! signal. 

ffl C . Monte Carlo Simulations 

A model of photon migration in tissues is necessary to define the relationship between the measured 
reflectance and the optical properties. Analytical solutions from the diffusion equation are not appropriate in 
our case because we are interested in the reflectance close to the source, at a distance comparable to the 
transport mean free path [mfp']. " We performed Monte Carlo simulations to predict the measured 
reflectance of an homogeneous semi-infmite turbid media. The code we used was extensively tested 2 ' • 22 . Any 
phase function can be implemented in discretized form. 
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Our simulations take into account the exact diameter of the illuminating and collecting fibers, as well as their 
numerical apertures (NA=0.28 in tissue). However, the distortion of the signal due to the size of the fibers 
we used (200 ran) has an almost negligible influence on the reflectance curve. 

The mismatch of index of refraction at the surface of the medium is also taken into account in our simulation, 
by using the Fresnel law for each photon reaching the surface. Simulations have been performed with the 
exact geometry of the probe as described in Fig.l, taking into account the mismatch of index of refraction 



^ D ' eWeeVmeVrb l b^ 

the air and the sample outside the probe (p > 1 .8 mm) (see Fig.3 and Fig.4). Nevertheless we found that the 
effect of the boundaries located at p<1.4 mm have a weak influence on the intensity collected by the fibers. 
To illustrate this, the reflectance obtained with the exact probe configuration is compared in Fig.4 to the 
ffj simplified case of the semi-infinite space (n^ = 1.4, n probe = 1.5). As expected, there are important 
P differences between the exact and simplified cases for p>1.5 mm, close to the limit of the probe. The 
Uj decreased reflectance at p > 1.8 for the exact case is due to the increase of internal reflection at the interface 
U1 between the medium and the air. However, the differences between the exact and simplified cases are less 
T than few percents for p < 1.4 mm, corresponding to the region where the measurements are performed. 
K The simplified case, assuming a cylindrical symmetry, is computationally much less time-consuming than the 
R "exact configuration". "Exact configuration" simulations require approximately ten times the number of 
fi photons to achieve statistical errors comparable for the simplified case. Following the result illustrated in 
Fig.4, we decided to employ only the simplified semi-infinite condition and restrict to p < 1.4 mm. 

3. RESULTS AND DISCUSSION 
A. Boundary effects 

During in vivo investigations the ideal case of a medium with a perfect plane surface is never realized. 
Moreover, the probe could be slightly pushed inside the tissue. Therefore the effect of the medium boundary 
on the measured reflectance is important to quantify experimentally. In Fig.5, we show two measurements in 
Intralipid (X. = 675 nm), one at the surface and the second inside the medium. The difference between these 
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two measurements is less than 5%, which could seem very surprising at first glance. However, as it was 
* shown by simulation in Fig.4, the boundary condition outside the probe has only a weak effect on the 
intensity measured by the six fibers. The boundary condition created by the probe itself, i.e. the index of 
refraction mismatch between the medium and the adhesive inside the probe, is much more critical. 
The negligible effect of the tissue boundaries shown in this section is an important advantage for clinical 
investigation. This experiment also clearly demonstrates that the sample volume investigated is principally 
confined to the region just beneath the probe surface. This point is developed further in the next section. 

B . Depth of tissue investigated 

To quantify more precisely the volume of tissue probed by our technique, the average depth of scattering 
~ events was recorded for each photon detected in Monte Carlo simulations (backscattered photons). To allow 

2 more general statements, distances are expressed in mfp' units. For typical turbid tissue 1 and near infrared 

3 wavelengths, jx/ is around 1 mm 1 and |i a less than 0.1 mm" 1 , which means that 1 mfp'~ 1mm. Fig.6 shows 
H3 that the average depth of scattering is approximately around 1 mfp*. Moreover it shows that for typical tissue 
^ optical properties, structures greater than 2 mfp', located beneath the probe are not likely to contribute 
H; significantly to the measured signal (for p < 1.5 mfp'). 

5 To evaluate experimentally the average depth probed, we performed experiments in Intralipid placing the 
d probe at the liquid surface, and moving an absorbing plate placed horizontally, as described in Fig.7. The 
" reflectance R(p, d) was measured for varying thickness d. The ratio R(p, d)/R(p, d=<») is reported in Fig.7. 
This figure shows that the intensity of the reflectance is decreased by approximately 20% if the medium is 2 
[mfp'] thick, and by 10% if it is 3 [mfp'] thick (for the albedo and the range of p considered). These results 
imply that 80% or 90% of the photons do not penetrate deeper than 2 [mfp'] or 3 [mfp'], respectively, into 
the medium. This confirms our simulation results (Fig.6) where the average depth of scattering was 
estimated to be approximately 1 mfp'. Thus, for typical biological tissues, our measurements are mainly 
sensitive to the region of tissue located within 2 mm of the surface and the investigated volume is on the 
order of 1 mm 3 . 
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C. Calibration and test on microsphere suspension. 

In order to perform absolute intensity measurements, calibration is performed on a solid turbid siloxane 
phantom of known optical properties (determined independently by frequency domain photon migration 2 3 ). 
As shown in Fig. 8, these measurements, multiplied by a factor independent of the fibers, fit well the Monte 
Carlo simulations performed with the phantom coefficients. This factor is defined as the calibration factor for 
a given wavelength. It is derived from me a surement s on the standard tissue phantom obtained at the end Q f_ 
each set of clinical measurements. The calibration was performed at the end of each set of clinical 
measurements. 

Experiments on microsphere suspensions (polystyrene sphere 0 1.072 ±0.019 Jim) were performed to 
assess the accuracy of our theoretical model and the calibration method. The scattering coefficient and the 
0 phase function of such turbid media can be precisely known using Mie theory 23 . As no dye was added to the 
IS suspension, ji a was considered to be equal to the water absorption. In Fig.8, a measurement of the 
f=j reflectance is compared to a simulation computed with the microsphere suspension coefficients (\i s 9 - 1 .0 
-J mm 1 , ^ a = 0.00041 mm" 1 , y = 2.2). The excellent agreement found here between experiments and simulation 
; confirms the accuracy of our simulations, and the validity of our calibration procedure. 

S D. Inverse problem 

m Our goal is to solve the inverse problem which consists of extracting optical coefficients from the reflectance 
data. The measurements of the reflectance intensity R(p) and the slope of InR(p), determined at a distance p 
= 1 mm, can be used to derive n 5 * and for a given y value. Fig.9 shows graphically the relationship 
between \i s ' and |l a and the two parameters R(p=l mm) and l3 p lnR(p=l mm)l. To illustrate the influence of 
the parameter y, two examples are superimposediy = 1.5 and y = 1.9. We see clearly in Fig.9 that \i* and \i a 
can not be determined uniquely from the two parameters R(p=l mm) and !9 p lnR(p=l mm) I if y is unknown. 
This indetermination may be resolved by the values of R(p) and/or l3 p lnR(p)l at other distances. Therefore 
the following procedure was defined for tissue measurements: 

(1) determination of \i g 9 and n a from R(p=l mm) and O p lnR(p=l mm)l for discrete values of y (for 
example y =1.0,1.5, 1.9, 2.2) 

ll 



(2) simulations with the different sets of u. s * and n a obtained 
■ (3) comparison between the simulations and the reflectance profile for distances 0.35<p<l .4 mm. 
This last step allows us to determine the value of y which gives the best fit. Points 1 to 3 can be done 
iteratively to evaluate y more precisely. The precision that can be obtained depends on the optical coefficients 
themselves, and on the experimental uncertainties. Nevertheless two important technical points should be 

noted here: — 
First, the determination of u/ is only weakly influenced by y, for u s ' close to 1 mm 1 . Indeed in Fig.9 the 

differences induced by y = 1.5 or y = 1.9 are typically 10% for . 

In contrast, absolute determination of u, is critically sensitive to by y. However, if y remains constant, 
relative variations of u a can be still precisely evaluated. This point is discussed with the results obtained of 
2 Intralipid measurements in the next section. Metabolism monitoring or drug monitoring could be therefore a 

H* potential application of such a probe. 

S Second, the experimental determination of l9 p lnR(p=l mm)l requires measurements at different distances 
| close to 1 mm. To minimize errors on l9 p lnR(p=l mm)l due to experimental artifacts, we perform a fit of the 
f 3 reflectance curve (0.5 mm <p 1.4 mm) with the function m lP "" exp(-m 3 p) which was always found to 
S fit well to Monte carlo simulations for this restricted range of distances (the same function was also proposed 
§ by Bolt and ten Bosh 24 ). The parameters R(p=l mm) and l3 p lnR(p=l mm)l are then derived from the fit. 
* Once the optical coefficients u a and u/ are derived, the validity of this procedure is double-checked, by 
comparing the curve obtained from the Monte Carlo simulation to the experimental profile. 

E. Phantom measurements 

We present in this section measurements on tissue-like phantoms (Intralipid™ and on microsphere 
suspensions). The value of y is a priori not known for Intralipid. Measurements of Intralipid phantoms with 
varying u, 4 and n 5 ' values were used to test the inversion procedure. Phantom optical properties were 
calibrated by the FDPM technique, (performed at large source-detector separations and therefore insensitive 
to y). For Intralipid measurements, reflectance profiles were found to be fit best when the parameter y was 
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between 1.6 and 1. 8 for k =674 to 849 nm. n a values derived from our method arc plotted in Fig. 10. a 
■ versus the values obtained by the FDPM technique. Fig. lO.a shows that 7=1 .8 leads to an overestimation of 
H a by approximately 0.005 mm' 1 , whereas 7=1.6 leads to an underestimation of |i a by approximately 0.01 
mm' 1 . This confirms that exact value of 7 is between 1.6 and 1.8. This example clearly illustrates that 
absolute values of |X a are sensitive to 7. However relative measurements of A|i a are achievable with a 
sensitivity of 0.005 mm' 1 using our system. The accuracy of the inversion procedure on n s ' is illustrated in 
Fig.lO.b. In the case of as already mentioned above, the influence df 7 is weaker than in the case of |i a . 
A difference of only 5% is found on the \x s J if 7=1.6 is used instead of 1.8. This is also approximately the 
variation in \i % y values obtained by multiple FDPM measurements. We also found that in the case of a 
-constant Intralipid concentration, the measured values of \i s y vary less than 2% when |! a is increased by 
q adding dye, which proves that the data inversion procedure effectively uncouples |i a and |A S \ 
Z Further assessments of the inversion procedure were performed on microsphere suspensions. In this case the 
Rvalue of 7 is known: 7=2.22. The optical coefficients derived from our local reflectance measurements are 
^compared to values obtained with a spatially-resolved transmittance method described elsewhere 2122 . 

Fig. 1 La and Fig. 11. b. show ^ a and m* values obtained by both methods: spatially resolved transmittance 
and reflectance. Phantoms made of microsphere suspensions and ink, at different concentrations, were used. 
An excellent correlation is found between the p. a and )!/ values, obtained by the two methods. Small 
systematic differences (typically 10%) are found when comparing absolute values. They are mainly due to 
errors occurring in the calibration procedures (for both methods). 

The relatively small errors we found here are typical of errors found when different techniques for measuring 
turbid media optical properties are compared 25 . Such errors could be avoided by multiple calibrations on 
several turbid samples of different known optical properties. However, the accuracy of tissue measurements 
is subject to other major limitations, due for example to their structure and heterogeneity 26 . Therefore, more 
accurate calibration is not necessarily required for absolue tissue measurements, since we are mainly 
interested in observing optical property differences between tissue physiological states. 
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F. In vivo measurements on brain tissues 

Clinical measurements of normal and malignant neural tissues were recorded in vivo during brain surgery 27 . 
Two different cases are reported here. Case 1 was a 3 year old male and case 2 was an 8 year old male. 
Different types of tissues were investigated in each case. Several measurements (typically 6) were always 
performed successively at a given location. The intensity fluctuations (typically on the order of 10%) for 
l ) i HStt measu r ements were mainly due to tissue heterogeneity an d sl i ght probe movements . The averag e_ 
reflectance was calculated for each location, as well as the standard deviation. Note that the uncertainty due to 
the apparatus, estimated from measurements on a phantom, are much lower (< 5%). Before each set of 
measurements, the blood from the surgical site was carefully irrigated away with saline, and the probe 
cleaned with a saline damped sponge. The measurements presented here were performed in parallel with 
□ frequency domain measurements (FDPM)" using a source-detector separation of 10 to 14 mm. 
5 Fig.l2.a and Fig.l2.b. show the measured parameters R(p=l mm) as a function of l3 p lnR(p=l mm)! 
m obtained for cases 1 and 2, respectively. These graphs are similar to Fig.9, except that the relationship 
jj; between the parameters R(p=l mm) and l9 p lnR(p=l mm)l and the optical coefficients p. a and H s ' is indicated 
: only qualitatively by two arrows for better clarity (different grids corresponding to different y values should 
3 be superimposed). Quantitative results are reported in Table 1 and Table 2. 

^ In case 1, measurements were performed on normal cerebral cortex (frontal lobe and temporal lobe), optic 
03 nerve astrocytoma (size- 1.3 cm) and normal optic nerve. In case 2, measurements were performed on the 
skull, deep cerebellar white matter with scar tissue (from a previous surgery), medulloblastoma (size~3.8 
cm), and deep cerebellar white matter (normal). Tumor dimensions were estimated from conventional 
imaging techniques (i.e. Computed Tomography and/or Magnetic Resonance Imaging). 
As discussed in section B, the depth probed is less than about 3 mm. For each type of tissue that we 
investigated, the influence of surrounding tissues on the measurement is weak. In particular, only gray matter 
is investigated during the cerebral cortical surface measurements. 

Fig.l2.a and Fig.l2.b. show that the parameters R(p=l mm) and l3 p lnR(p=l mm)l provide an excellent 
discrimination between tissue types. The spectroscopic signature on R(p=l mm) and l9 p lnR(p=l mm)l 
should also be noted. The cortex and the skull exhibit less significant spectroscopic differences compared to 
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tumor tissues such as the astrocytoma and the medulloblastoma. Thus, the parameters R(p=l mm) and 
* 0 p lnR(p=l mm)l could be useful for optical biopsy. Nevertheless we believe that in order to fully exploit 
these results, the differences found must be explained in terms of scattering and absorption parameters. 
These factors, in turn, can be used to understand physiological and structural variations. The procedure 
described in the section D was used to determine the coefficients u\ a and 7 from the measured curves. 
This procedure is fully illustrated for measurements of the temporal lobe and astrocytoma. Results of optical 

the FDPM technique are also indicated. 

Fig.l3.a shows that the best fit to the cortex (temporal lobe) data is obtained with 7=1 .9. Lower values, such 
as y=1.5, led to impossible values for R(p=l mm) and l9 p lnR(p=l mm)l m) and should be therefore 
m rejected. Larger values, such as 7=2.2 do not fit the reflectance data for distances p < 0.8 mm. The y value is 
£ then estimated to be y=1.9±0.2. The coefficients found with 7=1.9 are: u\ s ' = 1.0 mm' 1 and u\ a = 0.02 mm -1 at 
Q X = 674 nm, and u\/ = 0.82 mm 1 and |i a = 0.025 mm' 1 at X = 956 nm. Note that 7 = 2.2 would lead to 
jj| almost identical ji s ' (differences less than 5%) and an overestimation of |J. a of approximately 0.02 mm* . 
y Taking into account the uncertainties on the measurements (approximately 5% on R(p=l mm) and 
C ld p lnR(p=l mm)l for this tissue), and the uncertainty on 7 (±0.2), the error on the absolute value of |i s ' is 
H estimated to be 5%. 

;5 The astrocytoma reflectance obtained at 674 nm and 956 nm are plotted in Fig.l3.b. with simulations 

CD 

performed with 7 = 1.5 and 1.9. Both values of 7 can fit well the reflectance at distance p > 0.8. For shorter 
distances p < 0.8, the experimental curve falls between the simulations with 7=1.5 and 1 .9. The 7 value is 
then estimated to 1.7±0.2. Taking into account this uncertainty on 7, n s * = 1.25 ±0.10 mm' 1 , |X a = 0.14 
±0.03 nun 1 at X. = 811 nmand \i s 7 = 0.73 ±0.1 mm 1 , |i= 0.15±0.05 mm' 1 at X = 956 nm. As reported in 
Table 1, the absorption coefficients at 811 and 849 nm are lower than those obtained at X - 674 nm or >. 
=956 nm. This result is consistent with the fact that the main near-infrared tissue chromophores, hemoglobin 
and water, have absorption maximum at approximately at X < 700 nm and X = 970 nm, respectively 28 . The 
overall absorption is much higher in the tumor than the cortex presumably due to the greater hemoglobin 
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abundance of these components, particularly hemoglobin in the tumor. Tumors generally grow more new 
* blood vessels and thus have a higher blood flow. 

The |i s * of the tumor is similar to the cortex at X = 956 nm. However the variation of ji/ between X = 674 
nm and X - 956 nm is much larger for the tumor than for the cortex: p. 5 '(X=956 nm) - ^'(^=674 nm) « 
0.57 mm* 1 for tumor, n/(X=956 nm ) - |i s '(^=674 nm ) =0.10 mm" 1 for the cortex. Such spectroscopic 
variations may be attributable to structural differences between tissue types. Indeed such differences may 
depend on the average size or size distribution of scattering structures within or between cells. 
Higher scattering is found for the normal optic nerve compared to all the other tissues in case 1. This is likely 
due to the presence of myelin. In vitro measurements have also shown that myelin containing white matter 
exhibits a higher scattering coefficient compared to other tissue 1222 . The accuracy of the measurement on the 
m optic nerve tissue can be affected by the high anisotropy of this type of tissue. Indeed, it has been reported 
rr that light propagation depends on whether the direction considered is parallel or perpendicular to the nerve 
iTj fibers 29 " 31 . Interestingly, comparison between measurements with the two symmetric sources reveals larger 
{^differences than for the other tissues. This variation is represented in Fig.l2.a by the relatively large 
^ uncertainties associated with the optic nerve measurements. 

JZ In case 2, the first measurements were acquired directly from the skull. The difference in reduced scattering 
p; properties between X = 674 and 956 nm is very small (ji s * = 0.9 and 0.85 mm" 1 , respectively) as shown in 
jfiFig.l2.b. Previously determined in vitro values 32 for pig skull were approximately twice as large. Such 

We? 

differences may be due to variations in water content and sample preparation between in vivo and in vitro 
studies. 

Case 2 data were also obtained from, in series deep cerebellar white matter with scar tissue, tumor 
(medulloblastoma) and cerebellar white matter in the excavated tumor bed. Fig.l2.b. shows that scar tissue is 
well differentiated from other structures and is characterized by high y value (y = 2.2) and low absorption (^i a 
< 0.02 mm" 1 ), which is consistent with the low vascularization of such tissue. In contrast, larger |a a smaller y 
and large spectroscopic differences in are found in white matter and medulloblastoma. As in case 1, the 
large ji a values obtained for the tumor can be related to higher hemoglobin content, as is often found for 
cancerous tissues. Taking into account measurement variability, no significant differences are found between 
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the tumor and the "normal" white matter. However it is not clear that the so-called "normal" tissue, measured 

at the surgical boarder could be considered disease-free and unaffected by the tumor vasculature. 

The different type of tissue examined here showed and a clear dependence on the phase function parameter y. 

This confirms that y may be a valuable parameter for tissue characterization. Nevertheless, an improvement 

of the accuracy would be necessary to conclude about the importance of this parameter. 

The optical coefficients we obtained can be compared to the measurements performed simultaneously with 

"Mu^l^PM^ « 

compared to the spatially-resolved technique described here. Generally the |l s ' values are similar between the 
two methods, whereas more differences are found for |i a . For the cortex and skull measurements, both |i s * 
and |Li a values are in excellent agreement. The decrease of from X = 674 to 956 nm, is slightly more 
m pronounced in the FDPM data. In contrast, important differences are found for |a a results where tumor values 
P (astrocytoma and medulloblastoma) obtained by FDPM are significantly lower. This can be explained by the 
S sensitivity of the spatially-resolved probe to the high local hemoglobin content which can be resolved only by 
yi the small source-detector separations. In contrast, the large source detector separation employed by the 
~ FDPM probe interrogates much greater larger tissue volumes and hence measures average optical properties 
~ from multiple structures (e.g. normal + malignant). 

g 4. CONCLUSION 

The purpose of this work was to assess the performance of spatially resolved reflectance using short source- 
detector separations (< 1.4 mm). Monte Carlo simulations, the accuracy of which were confirmed by 
experiments on microsphere suspensions, were used to establish the correspondence between the measured 
reflectance and the optical properties. 

Optical properties which can be determined by this technique are the absorption coefficient |i a , the reduced 
scattering coefficient \i s \ and a parameter of the phase function rK^VO-g,), where g : (=g) and g 2 are the 
first and second moment of the phase function. Experiments on calibrated Intralipid solutions showed that n/ 
and \i & can be determined with a precision of ±0.05 mm" 1 and ±0.005 mm" 1 respectively, and y can be 
determined with a precision of typically 10%. Systematic errors are possible if the parameter y is not 
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determined with sufficient accuracy. These performances could be improved with by lowering the uncertainty 
* on reflectance measurements. 

Experiments and simulations helped to define the average volume probed by this technique. For typical 
tissues, the average probe depth is about typically 1 mm and the influence of layers located below 3 mm is 
negligible. 

Finally, in vivo measurements on human brain showed that excellent discrimination can be obtained between 
different types of neural tissues, normal and abnormal. Good correlation has been found between spatially- 
resolved reflectance and simultaneous measurements performed by frequency domain photon migration 
(FDPM). These two techniques offer interesting complementary features. The spatially resolved probe can 
potentially provide better differentiation between different types of tissue, due to its sensitivity to local 
m structure. This is due to the fact that substantially smaller volume of tissue is probed. On the other hand, due 
S to physical limitations imposed by large NIR mean absorption lengths in tissue, the precision for \i B estimate 
2 is likely to be worse. Consequently, the short distance, spatially-resolved technique appears to be well suited 
Hi for clinical settings, which require rapid localized tissue identification, such as endoscopic or needle-based 
^ "optical biopsy", and intraoperative tissue mapping for surgical guidance. 
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FIGURE CAPTIONS 

Fig. 1 . Probe for the measurement of the spatially resolved reflectance 
Experimental set-up. 

Boundary condition for exact probe configuration. 



Fig.2. 
Fig.3. 
Fig.4. 



Comparison between the exact and simplified boundary conditions (The numerical aperture was 
slightly increased to NA=0.37 instead of NA=0.28 for this graph, because of the small amount 
of detecied photons in the exact configuration) . The optical coefficient arc ; y . / — 1 mg 



Fig.5. 



1 Fig. 6. 



^Fig.7. 



=0.01 mm" 1 and 7=1.9. 

Comparison between measurements in a semi-infinite and an infinite media. The medium is 
Intralipid. The optical coefficients were measured by the FDPM technique \i g * = 1.2 mm" 1 , p. a = 

0.0005 mm" 1 (?t=675 nm). 

Average depth of scattering events was recorded for each photon detected in Monte Carlo 
simulations (backscattered photons). Case of y= 1.9. 

Effect of the Intralipid thickness on the reflectance. The set-up for the investigation on Intralipid 
with varying thickness is illustrated on right of the graph. The optical coefficients of the Intralipid 
were measured by the FDPM technique (k = 956 nm). The reduced albedo is a*=n s 7(m*+fi a ) = 

0.98 

Calibration measurement on a siloxane phantom and measurement test on a microsphere 
suspension. The calibration measurement is multiplied by the calibration factor to fit the 
corresponding simulation. The phantom optical properties were measured with frequency- 
domain photon migration: |X a =2.4± 0.2 10" 4 mm -1 , |X S '=1.82 ±0.007 mm" 1 (X=674 nm). The 

measurement on microsphere suspension is multiplied by the calibration factor derived from the 
phantom measurement, corresponding simulation is performed using the optical properties 

derived from Mie theory: |i/ = 1.0 mm" 1 , p. a = 0.0005 mm" 1 , Mie phase function (g=0 .916, y ~ 
2.2). 



! Fig.8. 



Fig.9. Relation between the parameters R(p=l mm) and l3 p lnR(p=l mm)l and the optical coefficients 
H,* and |x a Case of y =1.5 and 1.9. 
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Fig.lO.a. Comparison between n a obtained by the frequency domain photon migration technique (FDPM) 
- ' and by the probe measurements. Measurements on Intralipid and dye. 

Fig. lO.b. Comparison between \x* obtained by the frequency domain photon migration technique (FDPM) 
and by the probe measurements. Measurements on Intralipid 

Fig. 11. a. Comparison between ji a obtained from spatially-resolved transmittance (method described in 
references 21 and 22) and |i a from by the probe measurements. Measurements on microsphere 
suspension. 



Fig.ll.b. Comparison between JLI/ obtained from spatially-resolved transmittance (method described in 
references 21 and 22) and |i s * from by the probe measurements. Measurements on microsphere 
suspension. 

^Fig.l2.a. Clinical measurements in vivo on human brain. Case L R(p=l) and l3 p lnR(p=l mm)! for 

y '- 

O different types of brain tissues: normal cortex (frontal and temporal lobe), astrocytoma of optic 

^ nerve and normal optic nerve. 



HjFig.l2.b. Clinical measurements in vivo on human brain. Case 2. R(p=l) and l3_lnR(p=l mm)l for 

y] p 

^ different types of brain tissues: skull, deep cerebellar white matter, deep cerebellar white matter 

£^ with scar tissue and cerebellum tumor. 

ssa. 

I = 

i=j Fig. 13. a. Comparison between the spatially reflectance curve measured on normal cortex (temporal lobe) 
^ and simulations. 

Fig.l3.b. . Comparison between the spatially reflectance curve measured on astrocytoma of optic nerve and 
simulations. 
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TABLE CAPTIONS 

Table 1. Optical properties of normal and malignant human brain tissue. Casel 
Table 2. Optical properties of normal and malignant human brain tissue. Case2 
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Fig. 2 
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Fig. 4 
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Fig. 8 
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Fig. lO.b 
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Fig. 11-b 
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Fig. 12.a. 
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Fig. 12.b 
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